一 种 新 型 冠状 病毒 肺炎 (COVID-19) 药 物 与 天 然 产物 快速 发 现 的 计算 
药理 学 方法 


aCODE: Agile Discovery of Drugs and Natural Products for Emerging 
Epidemic COVID-19 based on Computational Pharmacology 


全 源 "?、 梁 峰 吉 *、 熊 江 交 E ?3 


"哈尔滨 工业 大 学 (深圳 ) 计 算 机 科学 与 技术 学 院 ， 广 东 深 圳 ，518055 

?深圳 市 绿 航 星际 太空 科技 研究 院 医 学 健康 部 表 观 遗传 学 实验 室 ， 广 东 深 圳 ，518117 
“航天 医学 基础 与 应 用 国家 重点 实验 室 ， 北 京 ，100094 

“Aromability Inc.Beijing,100080. 通讯 作者 :能 江 辉 , E-mail: xiongjh77@163.com 


[摘要 ] 新 发 传染 病 爆 发 流行 期 间 , 明 需 提出 候选 药物 功效 与 机 制 的 科学 假说 。 疫 
苗 或 新 药 研 发 均 需 要 一 定时 间 , 因而 药物 重 定位 ( 老 药 新 用 ) 策 略 有 其 独特 价值 。 
但 是 新 发 疾病 其 病原 体 、 宿 主 反 应 的 临床 数据 不 充分 , 制约 了 候选 药物 假设 的 提 
出 。 此 阶段 常 根 据 病人 临床 特征 进行 广 谱 抗 病毒 药物 的 尝试 。 本文 借鉴 人 工 智 能 
领域 常见 的 启发 式 搜索 思路 ， 提 出 一 种 新 方法 (aCODE)， 基 于 前 期 有 一 定 疗效 
提示 的 广 谱 抗 病毒 药 , 获得 其 宿主 靶 蛋 白 集 合 , 在 全 基因 组 尺度 上 搜索 与 之 相关 
性 最 高 的 基因 模块 组 合 ， 进 而 对 候选 化 合 物 (如 已 批准 上 市 药物 、 天 然 产 物 ) 进 
行 模式 匹配 与 统计 检验 排序 。 本 方法 可 根据 临床 实践 的 进展 更 新 输入 药物 , 迭代 
输出 更 精准 结果 , 输出 的 天 然 产物 或 中 药 、 药 食 同 源 成 分 结合 其 它 信 息 后 可 实施 
快速 测试 , 形成 敏捷 研发 测试 闭环 。 本 方法 的 第 二 版 更 新 及 其 与 文献 证 据 的 比 对 
分 析 请 参考 : http:/ /chinaxiv.org/abs/202002.00024。 
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Abstract 

During the outbreak of new infectious diseases, there is an urgent need to put 
forward scientific hypotheses for the efficacy，mechanism and side effects of 
candidate drugs. The research and development of vaccines or new drugs need a 
certain period of time, so the strategy of drug repositioning has its place. However, 
the clinical data of pathogen and host response of new diseases are not ready, 
restricts the hypothesis of candidate drugs. At this stage, we often try to use 
broad-spectrum antiviral drugs according to the clinical characteristics of patients. In 
this paper, we propose a new method aCODE (agile discovery method of drugs or 
natural products for emerging epidemic) which based on the heuristic search strategy 
widely used in the field of artificial intelligence. Based on the broad-spectrum 
antiviral drugs with some early efficacy tips, the host target protein collection is 
obtained, and the associated gene modules is searched on the whole genome scale. 
We then carry out pattern matching and statistics for candidate compounds (such as 


1 


approved drugs and natural products ingredients). This method can update the input 
drugs according to the progress of clinical practice, and output more accurate results 
iteratively. The output components from natural products, traditional Chinese 
medicine or food can be used for quick trial to form a closed loop of agile R & D test. 
In addition, for the second update of this method and its comparison with literature 
evidence, please refer to: http://chinaxiv.org/abs/202002.00024. 

Key words: COVID-19; computational pharmacology; drug repositioning; natural 
products; Network pharmacology 
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重大 新 发 传染 病 爆 发 期 的 特点 是 时 间 紧 迫 、 数 据 有 限 。 基 于 病毒 靶 标 的 药物 
发 现 路 径 需 要 对 病毒 结构 与 功能 的 精细 了 解 , 而 临床 疗效 毒 副作用 验证 需要 大 量 
时 间 。 基 于 表 型 的 药物 重 定位 ( 老 药 新 用 ) 技术 是 一 种 重要 的 新 思路 [1]， 但 是 
对 于 新 发 传染 病 而 言 ， 短 期 很 难 获得 足 量 的 基因 组 学 、 表 型 组 学 数据 。 针 对 这 种 
数据 极其 有 限 情 况 下 需要 快速 产生 可 验证 药物 假说 并 进行 量化 评价 的 需求 , 我 们 
崩 鉴 人 工 智能 领域 常用 的 启发 式 搜索 策略 ， 提 出 一 种 新 方法 即 aCODE (agile 
discovery method of drugs or natural products for emerging epidemic) 
方法 (图 1) 。 基 于 有 限 疗 效 的 广 谱 抗 病毒 药 ， 获 得 其 宿主 轰 蛋 白 集 合 ， 在 全 基 
因 组 尺度 上 搜索 与 之 相关 性 最 高 的 基因 模块 组 合 , 进而 对 候选 化 合 物 (如 已 批准 
上 市 药物 、 天 然 产物 ) 进行 模式 匹配 与 统计 检验 排序 。 本 文 以 最 近 在 我 国 爆发 的 
新 型 冠状 病毒 COVID-19 为 例 ， 以 部 分 报道 有 一 定 疗效 提示 的 药物 为 输入 ， 对 
天 然 产物 成 分 进行 计算 排序 , 并 将 初始 药物 与 天 然 产 物 分 子 混合 聚 类 , 以 了 解 候 
选 化 合 物 发 挥 作用 的 潜在 细胞 分 子 机 制 。 本 方法 可 以 随 临 床 实 践 的 进展 将 初始 药 
物 列表 进行 迭代 更 新 ， 从 而 优化 系统 的 精准 度 与 预测 能 力 ; 输出 的 天 然 产 物 或 中 
药 、 药 食 同 源 成 分 可 实施 快速 测试 ， 形成 敏捷 研发 测试 闭环 ; 并 为 爆发 期 医疗 资 
源 有 限 情况 下 大 量 人 群 的 居家 预防 与 方案 测试 提供 科学 支持 。 
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图 1. aCODE (agile discovery method of drugs or natural products for 
emerging epidemic) 方 法 的 框架 及 主要 分 析 流 程 


2. 数据 资源 与 分 析 方 法 


2.1 人 类 基因 模块 的 鉴定 

在 本 研究 中 ， 人 类 基因 模块 数据 源 下 载 自 Molecular Signatures Database 

(MSigDB) 数据 库 (https://www.gsea-msigdb.org/gsea/msigdb/) [2]。 数 据 库 包含 

常用 的 GO BiologicalProcess，GO CellularComponent，GO MolecularFunction 等 
基因 模块 。 

上 述 基因 模块 信息 是 静态 的 ,是 参与 某 生物 学 过 程 的 所 有 基因 的 集合 
衰老 是 影响 人 类 各 类 疾病 、 免 疫 功 能 变化 的 最 普遍 因素 ,筛选 与 年 龄 ， 
作为 各 个 模块 的 代表 。 采 用 包含 年 龄 信息 的 人 类 血细胞 DNA 甲 基 化 数据 
GSE40279 (https:/www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE40279) [3]， 
该 数据 集 包含 来 自 人 类 全 血 的 656 个 样本 的 DNA 甲 基 化 谱 , 年 龄 从 19 岁 到 101 
岁 不 等 。 利 用 新 近 提 出 的 可 提取 基因 水 平 DNA 甲 基 化 特征 的 算法 [4-6]， 筛 选 衰 
老 相 关 的 、 受 DNA 甲 基 化 调控 的 基因 集合 。 该 算法 (SIMPO) 以 每 个 基因 的 
SIMPO 得 分 来 代表 基因 水 平 的 DNA 甲 基 化 程度 。 

SIMPO 算法 的 计算 公式 如 下 : 


SimPo score = 


ym+n— 2) (1-1) 
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其 中 ， 


S， 三 rr —1)Sf + (n— 1)S? 
:位 于 基因 body 区 域 的 所 有 探 针 的 平均 DNA 甲 基 化 值 ; 
位 于 基因 promoter 区 域 的 所 有 探 针 的 平均 DNA 甲 基 化 值 ; 
: 位 于 基因 body 区 域 的 所 有 探 针 数目 ; 
位 于 基因 promoter 区 域 的 所 有 探 针 数目 ; 
3: 位 于 基因 body 区 域 探 针 的 DNA 甲 基 化 值 的 方差; 
52: 位 于 基因 promoter 区 域 探 针 的 DNA 甲 基 化 值 的 方差 。 
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接着 ， 本 研究 基于 工 检验 计算 了 老年 (>=80 岁 ) 和 青年 (<35 岁 ) 个 体 之 间 
的 基因 SIMPO 得 分 差异 , 挑选 工 检验 结果 中 了 值 <“0.001 的 基因 作为 年 龄 显著 相 
关 基 因 模 块 。 该 模块 与 MSigDB 数据 库 收集 的 各 个 生物 学 功能 基因 模块 的 交集 
将 作为 各 个 生物 学 功能 基因 模块 的 代表 。 由 此 , 本 研究 获得 了 多 个 与 衰老 相关 的 
生物 学 功能 基因 模块 。 
2.2 初始 药物 显著 关联 基因 模块 鉴定 

当前 版 本 ， 输 入 表 1 的 10 个 药物 作为 初始 药物 ， 主 要 来 源 于 国家 第 四 版 指 
南 和 部 分 有 细胞 学 证 据 的 药物 [7]， 及 广 谱 抗 病 毒药 。 


表 1. 初始 药物 集合 


Drug name 药物 名 称 Drug name 药物 名 称 
Lopinavir 罗 匹 那 韦 Ribavirin 利 巴 韦 林 
Ritonavir 利 托 那 韦 Indinavir 节 地 那 韦 
Oseltamivir 奥 司 他 韦 Saquinavir 沙 奈 那 韦 
Chloroquine 握 众 Carfilzomib 卡 非 佐 米 
Mycophenolate mofetil 霉 酚 酸 酯 Nelfinavir 奈 非 那 书 


上 述 抗 病毒 初始 药物 的 邯 标 收集 自 STITCH (chemical association networks， 
http://stitch.embl.de) 数据 库 [8]。 为 了 兼顾 STITCH 数据 库 中 靶 标 数量 较 少 的 化 
合 物 ， 本 研究 只 选择 每 个 化 合 物 综合 打分 (combined score) 前 300 的 靶 标 作为 
该 化 合 物 的 候选 靶 标 。 基 于 超 几 何 分 布 检验 来 计算 初始 药物 与 各 个 功能 基因 模块 
的 相关 性 。 挑 选 超 几何 分 布 检 验 结果 中 P 值 最 显著 top200 的 生物 学 功能 基因 模 
块 作为 初始 药物 显著 关联 基因 模块 。 

超 几 何 分 布 检验 公式 如 下 : 


P—Value=1—F(x—1/M,K,N)=1— 7 2) 
其 中 x 是 初始 药物 与 功能 基因 模块 交集 的 基因 数量 ，M 是 全 集 基因 数量 ，K 
是 初始 药物 靶 标 的 基因 数量 ，N 是 某 个 生物 学 功能 基因 模块 的 基因 数量 。 


2.3 潜在 抗 新 型 冠状 病毒 COVID-19 药物 预测 
本 研究 用 于 筛选 抗 新 型 冠状 病毒 药物 的 背景 药物 数据 集 收 集 自 天 然 产 物 (中 

药 ) 化 学 成 分 数据 库 Traditional Chinese Medicine Integrated Database (TCMID， 
http://www.megabionet.org/tcmid/) [8]。 与 抗 新 型 冠状 病毒 COVID-19 初始 药物 一 
致 ， 背 景 药物 的 靶 标 收集 自 STITCH 数据 库 ， 并 选择 每 个 化 合 物 综合 打分 

(combined score) 前 300 的 诅 标 作为 该 化 合 物 的 候选 靶 标 。 接 着 ， 通 过 超 几 何 
分 布 检验 , 我 们 分 别 计算 了 背景 药物 与 各 个 初始 药物 显著 关联 基因 模块 的 相关 性 。 
此 外 , 基于 初始 药物 与 背景 药物 计算 的 超 几 何 分 布 检验 P 值 , 本 研究 通过 聚 类 方 
法 获得 了 初始 药物 与 背景 药物 的 相关 性 。 


3. 结果 


本 研究 基于 SIMPO 算法 , 获得 了 年 龄 相关 的 DNA 甲 基 化 数据 中 每 个 基因 的 
SIMPO 得 分 来 代表 基因 水 平 的 DNA 甲 基 化 程度 。 基 于 工 检验 , 本 研究 计算 了 老 
年 (>=80 岁 ) 和 青年 (<35 岁 ) 个 体 之 间 的 基因 SIMPO 得 分 差异 ， 挑 选 工 检验 
结果 中 P 值 <0.001 的 基因 作为 年 龄 显著 相关 基因 模块 。 接 着 ， 通 过 MSigDB 数 
据 库 中 生物 学 功能 基因 模块 与 年 龄 显著 相关 基因 模块 做 交集 , 本 研究 鉴定 了 衰老 
相关 的 生物 学 功能 基因 模块 。 

接着 , 结合 抗 病 毒 初始 药物 的 靶 标 基因 及 衰老 相关 的 生物 学 功能 基因 模块 数 
据 , 通过 超 几何 分 布 检验 获得 P 值 ,本 研究 鉴定 了 200 个 与 抗 病毒 初始 药物 显著 
关联 的 生物 学 功能 基因 模块 (图 2)。 其 中 ,图 3 显示 的 是 抗 病毒 初始 药物 Ribavirin 
显著 关联 的 生物 学 功能 基因 模块 。 
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图 2. 与 COVID-19 具有 显著 相关 性 的 基因 模块 (Top15) 
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图 3. 抗 病毒 初始 药物 Ribavirin 显著 关联 的 生物 学 功能 基因 模块 (Top15) 


接着 ， 本 研究 对 天 然 产物 (中药) 化 学 成 分 进行 抗 病毒 药物 筛选 。 结 合 天 然 
产物 〈 中 药 ) 化 学 成 分 在 STITCH 数据 库 中 的 靶 标 基因 ， 通 过 超 几 何 分 布 检验 ， 
我 们 分 别 计算 了 候选 药物 与 各 个 初始 药物 显著 关联 基因 模块 的 相关 性 得 分 , 当 一 
个 化 合 物 对 应 多 个 初始 药物 显著 关联 基因 模块 , 则 对 -log(P 值 ) 进 行 累加 ( 附 表 1)， 


其 中 p 值 是 超 几何 分 布 检验 的 p 值 。 


分 析 结果 显示 : 初始 输入 药物 Ribavirin 与 COVID-19 疾病 相关 基因 模块 的 关 
联 性 最 强 〈 附 表 1) 。 此 外 ， 有 4 个 天 然 产物 成 分 的 相关 性 得 分 高 于 已 知 的 抗 病 
毒 初始 药物 (包括 二 氧化 硅 ，3,4- 茶 并 荡 ， 芹 黄 素 ， 富 马 酸 ) 。 值 得 注意 的 是 ， 
相关 性 得 分 排名 靠 前 的 芦 丁 (rutin)、 芹 菜 素 (apigenin) 与 山 泰 酚 (Kaempferol) ， 
在 国家 卫 健 委 及 各 中 医 机 构 推 荐 中 药 处 方 草药 中 出 现 频次 较 高 ( 表 2) ,提示 本 
方法 有 利于 促进 潜在 抗 新 型 冠状 病毒 COVID-19 药物 的 发 现 。 


此 外 ， 基 于 初始 药物 与 候选 药物 计算 的 超 几何 分 布 检 验 P 值 ， 通 过 药物 - 基 
因 模 块 双 向 聚 类 , 可 以 观察 候选 药物 与 初始 药物 在 基因 模块 相关 性 上 的 相似 程度 
(图 4 和 图 5) 。 与 已 知 的 初始 药物 聚 类 较 近 的 候选 天 然 产物 化 合 物 将 被 预测 为 
潜在 的 抗 新 型 冠状 病毒 COVID-19 药物 ， 并 且 具 有 类 似 的 作用 机 制 (Mechanisms 


of Action)。 


图 4 提示 候选 药物 与 初始 药物 作用 机 制 的 育 关 图 《药物 视角 ) 


图 5. 提示 候选 药物 与 初始 药物 作用 机 制 的 聚 类 图 (基因 模块 视角 ) 


4 讨论 


本 文 结果 提示 ， 基 于 新 发 传染 病 早 期 有 限 的 用 药 经 验 反馈 ， 推 测 并 计算 疾病 
特征 进而 对 更 大 范围 的 候选 药物 进行 量化 分 析 是 可 行 的 。 本 方法 的 实质 是 先 在 药 
物 这 个 维度 空间 , 聚集 一 组 有 一 定 疗 效 提示 的 药物 ; 用 此 信息 指导 在 基因 维度 空 
间 上 计算 各 个 基因 模块 与 该 疾病 的 相关 性 并 进行 基因 模块 的 聚 类 ; 最 后 回 到 药物 
空间 对 候选 化 合 物 进 行 聚 类 与 排序 。 


表 2. 国 家 卫 健 委 及 各 中 医 机 构 推荐 中 药 处 方 中 草 药 化 学 成 分 出 现 频次 表 


化 合 物 
Linalool 
Stigmasterol 
Vitamin C 
Gamma-Sitosterol 
Camphor 
Eugenol 
Kaempferol 
Rutin 
Palmitic Acid 
Quercetin 
Adenine 
Camphene 
Alpha-Pinene 
Daucosterol 
Gallicacid 
Limonene 
Stearic Acid 


频次 
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化 合 物 
Sucrose 
D-Limonene 
Geraniol 
Naphthalene 
Vitamin B 
Apigenin 
Hexadecane 
Cholesterol 
Citral 
Higenamine 
Emodin 
Acetic Acid 
Scopoletin 
Thymol 
Stearicacid 
Phytol 
Vitamin B1 


频次 


32 
31 
31 
31 
31 
30 
30 
29 
29 
29 
28 
27 
27 
27 
26 
25 
29 


化 合 物 
Anethole 
Beta-Elemene 
Lauric Acid 
Trigonelline 
Alpha-Limonene 
Benzaldehyde 
Catffeicacid 
Campesterol 
Oleic Acid 
Berberine 
Beta-Sitosterol 
Choline 
Chrysophanol 
Dibutyl Phthalate 
Ecdysterone 
Eugenol Methyl Ether 


频次 
24 
24 
24 
24 
23 
23 
23 
23 
23 
22 
22 
22 
22 
22 
22 
22 


本 文 提出 的 aCODE 方法 具有 如 下 特点 : (1) 敏捷 迭代 性 : 本 方法 可 以 随 临 
床 实践 的 进展 将 初始 药物 列表 进行 迄 代 更 新 , 从 而 优化 系统 的 精准 度 与 预测 能 力 ; 
(2) 快速 闭环 性 : 输出 的 天 然 产 物 或 中 药 、 药 食 同 源 成 分 可 实施 快速 测试 ， 形 
成 敏捷 研发 测试 闭环 ; (3) 药理 机 制 方面 ， 可 对 候选 药物 与 初始 药物 的 药物 作 
用 的 细胞 分 子 生物 学 机 制 进 行 聚 类 、 比 对 。 并 为 爆发 期 医疗 资源 有 限 情况 下 大 量 
人 群 的 居家 预防 与 方案 测试 提供 科学 支持 。 

aCODE 方法 本 身 可 以 持续 更 新 的 方面 包括 : (1) 药物 靶 标 数据 库 的 更 新 。 
本 研究 目前 使 用 的 药物 部 标 数 据 来 自 公开 的 STITCH 数据 库 , 后 续 可 以 接 入 更 精 
细 的 机 构 专 有 药物 彼 标 数据 库 ; (2) 疾病 相关 基因 的 更 新 。 如 果 有 病毒 感染 的 
宿主 细胞 组 学 等 数据 , 可 以 对 本 文 用 到 的 疾病 相关 基因 模块 数据 进行 更 新 ; (3) 
基因 模块 特异 性 的 更 新 。 当 前 所 用 的 MsigDB 基因 模块 ， 大 部 分 来 自 Gene 
Ontology， 该 知识 表示 体系 是 一 个 层级 体系 ， 为 了 兼顾 特异 性 与 一 般 性 ， 选 择 了 
基因 数量 在 50-500 个 之 间 的 基因 模块 ， 后 续 可 以 对 该 数据 库 进行 扩充 或 微调 ， 
并 且 优 先 选择 基因 数量 比较 小 、 特 异性 更 高 的 基因 模块 ; (4) 其 他 药物 证 据 的 
整合 能 力 。 本 方法 是 一 个 纯 数据 驱动 的 计算 过 程 , 如 果 该 结果 与 其 他 来 源 的 临床 
或 文献 证 据 相符 ， 按 照 贝 叶 斯 原理 ， 则 我 们 对 该 候选 物 的 信念 与 兴趣 增强 ， 后 续 
可 以 结合 其 他 因素 进一步 测试 其 成 药 潜力 。 

值得 注意 的 是 ， 本 方法 是 一 个 类 似 GPS 的 快速 聚焦 锁定 候选 化 合 物 的 工具 ， 
具体 到 每 个 候选 化 合 物 , 需要 综合 各 个 维度 考察 其 发 挥 作用 的 大 小 、 方 向 、 毒 前 
作用 、 成 药性 、 经 济 性 等 。 本 方法 只 涉及 药物 与 候选 化 合 物 的 计算 药理 学 分 析 ， 
不 对 具体 的 候选 药物 进行 临床 提示 。 
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附 表 1。 COVID-19 药物 与 天 然 产 物 来 源 的 候选 化 合 物 及 其 分 子 机 制 


Compound 名 称 类 相关 显著 关联 基因 模块 (Top 3) 
name 型 性 
Ribavirin 利 巴 韦 林 D 377 GSE13485_PRE_VS_POST_YF17D_VACCINATION_PBMC_DN, 


GSE19888_ADENOSINE_A3R_INH_VS_ACT_WITH_INHIBITOR_PRETREATMENT 


_IN_MAST_CELL_UP GSE42724_NAIVE_ BCELL VS_ PLASMABLAST_UP 


Silicon 二 氧化 硅 327 GO_COFACTOR_METABOLIC PROCESS, GO_CYTOKINE_ACTIVITY, 
Dioxide GO_REGULATION_OF_IMMUNE_ EFFECTOR_PROCESS 
3,4-Benzop 3,4- 茶 并 251 GO_DRUG_ METABOLIC PROCESS, 
yrene 花 GO_OXIDOREDUCTASE ACTIVITY_ACTING_ON_PAIRED_DONORS_WITH_INC 


ORPORATION_OR_REDUCTION_OF_ MOLECULAR_OXYGEN, 
GO_COFACTOR_METABOLIC PROCESS 
Apigenin 芹菜 素 247 GO_MAMMARY_GLAND_DEVELOPMENT, PID_P53_DOWNSTREAM._ PATHWAY, 
GO_ANATOMICAL_ STRUCTURE HOMEOSTASIS 
Chloroquine 氛 唑 D 225 GO_DRUG METABOLIC_ PROCESS, GO_COFACTOR_METABOLIC_PROCESS, 
GO_PRODUCTION_OF MOLECULAR_ MEDIATOR_OF_ IMMUNE_ RESPONSE 
Fumaricacid 富 马 酸 216 GO_PROTEIN_TYROSINE_ KINASE ACTIVITY， 
GO_PROTEIN_AUTOPHOSPHORYLATION, 
REACTOME_SIGNALING_BY_RECEPTOR_TYROSINE KINASES 


Mycopheno ”和 霉 酚 酸 酯 D 204 GO_ENDOPEPTIDASE_ACTIVITY, GO_SERINE_HYDROLASE_ACTIVITY, 
late Mofetil GO_HUMORAL IMMUNE_RESPONSE 
1-Methoxyc ”1- 甲 氧 鸟 204 MODULE_197, KEGG_ERBB_SIGNALING_ PATHWAY, 
anthin-6-On ” 氨 酸 -6 一 PID_P53_DOWNSTREAM_PATHWAY 
e 
Atenolol 阿 蔡 洛 尔 197 GO_PROTEIN_TYROSINE_KINASE_ACTIVITY, 


GO_PROTEIN_AUTOPHOSPHORYLATION, 
GO_POSITIVE_REGULATION_OF_ PHOSPHATIDYLINOSITOL, 3_KINASE SIGNA 


LING 
L-Homocys L- 同 型 半 194 GO_DRUG METABOLIC_ PROCESS, GO_COFACTOR_METABOLIC_PROCESS, 
teine 胱 氨 酸 GO_COLLAGEN_CONTAINING_EXTRACELLULAR_MATRIX 
L-AlphaGa L-a, Y- 190 GO_PROTEIN_TYROSINE_ KINASE ACTIVITY， 
mma-Diami 二 氨基 ] GO_PROTEIN_AUTOPHOSPHORYLATION, 
nobutyric 酸 GO_POSITIVE_REGULATION_OF_ PHOSPHATIDYLINOSITOL 3_KINASE SIGNA 
Acid LING 
Baicalein 黄 芬 黄 素 190 GO_PEPTIDYL_THREONINE_MODIFICATIION， 


KEGG_ERBB_SIGNALING_ PATHWAY, 
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Compound 名 称 类 相关 显著 关联 基因 模块 (Top 3) 
name 型 性 


GO_BLOOD_VESSEL_ ENDOTHELIAL_CELL_ MIGRATION 


1,7-Diphen ”1,7- 二 茶 188 REACTOME_NEGATIVE_REGULATION_OF_THE,_ PI3K_AKT_NETWORK, 
yl-5-Hepten ” 基 -5- 庚 烯 GO _ MAMMARY_GLAND_DEVELOPMENT, 
e-3-One 3 GSE7509_UNSTIM_VS_FCGRIIB_STIM_DC_DN 
Benzyl 水 杨 酸 他 188 REACTOME NEGATIVE_REGULATION_OF_THE_PI3K_AKT_NETWORK， 
Salicylate 酯 GO_MAMMARY_GLAND_ DEVELOPMENT, 


GSE7509_UNSTIM_VS_FCGRIIB_STIM_DC_DN 
Creatinine 肌 酝 187 MODULE_24, GO_CYTOKINE_ACTIVITY, 
GO_REGULATION_OF_ VASCULATURE DEVELOPMENT 


Benzene 葵 186 GO_DRUG_METABOLIC_ PROCESS, PID_P53_DOWNSTREAM_PATHWAY, 
GO_OXIDOREDUCTASE_ACTIVITY_ACTING_ON_PAIRED_DONORS_WITH INC 
ORPORATION_OR_REDUCTION_OF MOLECULAR_OXYGEN 
Palmitoleica ”棕榈 油 酸 179 GO_PROTEIN_TYROSINE_KINASE_ACTIVITY, 
cid GO_PROTEIN_AUTOPHOSPHORYLATION, 
GO_OXIDOREDUCTASE_ACTIVITY_ACTING_ON_PAIRED DONORS_WITH_INC 
ORPORATION_OR_REDUCTION_OF MOLECULAR_OXYGEN 
2-Propenal 。 2- 丙 烯 醛 169 GO_ENDOPEPTIDASE_ACTIVITY, GO_SERINE_HYDROLASE_ACTIVITY 
KEGG_ERBB_SIGNALING_PATHWAY 
Neoisorutin ”新 异 芦 丁 166 GO_ENDOPEPTIDASE_ACTIVITY, GO_SERINE_HYDROLASE_ACTIVITY, 
REACTOME_NEGATIVE_REGULATION_OF_THE,_ PI3K_AKT_NETWORK 
Rutin 芦 丁 166 GO_ENDOPEPTIDASE_ACTIVITY, GO_SERINE_HYDROLASE_ACTIVITY, 
REACTOME_NEGATIVE_REGULATION_OF_THE,_ PI3K_AKT_NETWORK 
Deoxycholi ”脱氧 胆 酸 165 GO_OXIDOREDUCTASE_ACTIVITY_ACTING _ON_PAIRED DONORS_WITH_INC 
cacid ORPORATION_OR_ REDUCTION_OF MOLECULAR_OXYGEN, 


GO_DRUG METABOLIC_ PROCESS, KEGG_ ERBB_SIGNALING PATHWAY 


[evil 


Delphinidin ”德尔 菲 157 REACTOME_NEGATIVE_REGULATION_OF_THE,_PI3K_AKT_NETWORK, 
丁 GO_REGULATION_OF PEPTIDASE_ACTIVITY， 
GSE32986_CURDLAN_HIGHDOSE_VS_GMCSF_AND_CURDLAN_HIGHDOSE_STI 
M_DC_DN 
Zerumbone 花 姜 酮 157 GO_ENDOCRINE_SYSTEM_DEVELOPMENT, 
GO_PEPTIDYL_THREONINE_MODIFICATION, 


ZHOU_INFLAMMATORY_RESPONSE_LPS_UP 


Wedelolacto ” 韦 地 内 酶 157 NABA_SECRETED_FACTORS, GO_CYTOKINE_ACTIVITY, 
ne GO_ENDOCRINE_SYSTEM_DEVELOPMENT 
Anwulignan ” 安 五 脂 素 157 GO_PHAGOCYTOSIS, GO_CYTOKINE_ACTIVITY, 
RUAN_RESPONSE_ TO_TNF_DN 
Macelignan ” 肉 豆 获 木 157 GO_PHAGOCYTOSIS, GO_CYTOKINE_ACTIVITY, 
酚 素 RUAN_RESPONSE_TO_TNF_DN 
Ascorbic 抗坏血酸 152 GO_COFACTOR_METABOLIC_PROCESS， 
Acid GO_OXIDOREDUCTASE_ACTIVITY_ACTING_ON_PAIRED DONORS_WITH INC 


ORPORATION_OR_REDUCTION_OF_ MOLECULAR_OXYGEN, 
MOOTHA_MITOCHONDRIA 
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Compound 
name 


Threonine 


Methylpred 


nisolone 


Kaempferol 


Anthracene 


Piceatannol 


Eckol 


3,4,5-Trihyd 
roxybenzoic 
Acid 


Gallicacid 


1-Methoxy- 
4-(1-Propen 
yl)-Benzene 
1-Methoxy- 
4-(1-Propen 
yDBenzene 


Anethole 


Cis-Anethol 


e 


P-Propenyl- 


Anisole 


Trans-Aneth 


ole 


名 称 


山 泰 栈 


埃 科 尔 


3,4,5- 三 羟 


会 


基 茶 甲酸 


没食子 酸 


1- 甲 氧 基 
-4-(1- 丙 烯 


基 )- 莱 


1- 甲 氧 基 
-4-(1- 丙 烯 
基 ) 葵 
茄 香 脑 


顺 式 苘 香 
脑 


对 丙 烦 基 


A 
[三 
薄 


类 ”相关 
型 性 


150 


150 


148 


147 


144 


144 


143 


143 


139 


139 


139 


139 


139 


139 


显著 关联 基因 模块 (Top 3) 


GO_ANTIGEN_RECEPTOR_MEDIATED_SIGNALING PATHWAY, 
GO_IMMUNE_ RESPONSE_ REGULATING CELL_ SURFACE RECEPTOR_SIGNALI 
NG_PATHWAY, GO_PEPTIDYL_THREONINE MODIFICATION 
GO_ADAPTIVE_IMMUNE_RESPONSE_BASED_ ON_SOMATIC RECOMBINATION 
_OF_IMMUNE_ RECEPTORS_BUILT_FROM_ IMMUNOGLOBULIN_SUPERFAMILY 
_DOMAINS, GO_LYMPHOCYTE_MEDIATED IMMUNITY, 
SMID_BREAST _ CANCER NORMAL LIKE UP 
GO_OXIDOREDUCTASE_ACTIVITY_ACTING_ ON_PAIRED DONORS_WITH_INC 
ORPORATION_OR_REDUCTION_OF_ MOLECULAR_OXYGEN, 

GO_ MAMMARY_GLAND _ DEVELOPMENT, 
GO_BLOOD_VESSEL_ ENDOTHELIAL_CELL _ MIGRATION 
GO_CATALYTIC_ACTIVITY_ACTING_ON_RNA, 
GO_CAIALYTIC_ACTIVITY_ACTING_ON_DNA， 

GO_DNA CONFORMATION_CHANGE 
GO_PROTEIN_TYROSINE_ KINASE ACTIVITY， 
GO_REGULATION_OF_PEPTIDASE_ACTIVITY, 
GO_NEGATIVE REGULATION_ OF HYDROLASE_ ACTIVITY 
GO_ENDOCRINE_SYSTEM_DEVELOPMENT, 
KEGG_ERBB_SIGNALING_ PATHWAY, 
GO_PEPTIDYL_THREONINE MODIFICATION 
GO_REGULATION_OF_PEPTIDASE _ACTIVITY, 
PID_P53_DOWNSTREAM._ PATHWAY, GO_MAMMARY_GLAND DEVELOPMENT 


GO_REGULATION_OF_PEPTIDASE_ACTIVITY, 
PID_P53_DOWNSTREAM._ PATHWAY, GO _MAMMARY_GLAND_ DEVELOPMENT 
GO_COFACTOR_METABOLIC_PROCESS, 
GO_ENDOCRINE_SYSTEM_DEVELOPMENT, 
KEGG_ERBB_SIGNALING PATHWAY 
GO_COFACTOR_METABOLIC_PROCESS， 
GO_ENDOCRINE_SYSTEM_DEVELOPMENT, 
KEGG_ERBB_SIGNALING PATHWAY 
GO_COFACTOR_METABOLIC_PROCESS, 
GO_ENDOCRINE_SYSTEM_DEVELOPMENT, 
KEGG_ERBB_SIGNALING PATHWAY 
GO_COFACTOR_METABOLIC_PROCESS， 
GO_ENDOCRINE_SYSTEM_DEVELOPMENT, 
KEGG_ ERBB_SIGNALING PATHWAY 
GO_COFACTOR_METABOLIC_PROCESS, 
GO_ENDOCRINE_SYSTEM_DEVELOPMENT, 
KEGG_ERBB_SIGNALING_ PATHWAY 
GO_COFACTOR_METABOLIC_PROCESS， 
GO_ENDOCRINE_ SYSTEM_DEVELOPMENT, 
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Compound 名 称 
name 

3354.57 -33 

-Hexahydro 7- 六 羟基 
xyflavone 黄酮 
Myricetin 杨梅 素 
Honokiol 厚 朴 酚 

Dihydrotest ”二 氧 摆 酮 

OSterone 

4-Methylsul ”4- 甲 基 亚 
finyl Butyl 磺 酰 基 
Isothiocyan ” 基 异 硫 氰 

ate 酸 酯 
Sulforathan 硫 氟 烷 
ee 
注 ， “类型” 栏 标 当 


类 ”相关 


型 ”性 


137 


137 


134 


133 


133 


133 


显著 关联 基因 模块 (Top 3) 


KEGG_ ERBB_SIGNALING PATHWAY 
GO_POSITIVE_REGULAIION_OF_PHOSPHATIDYLINOSITOL 3_KINASE SIGNA 
LING, KEGG_ERBB_SIGNALING_ PATHWAY, 
GO_ENDOTHELIAL_CELL PROLIFERATION 
GO_POSITIVE_REGULATION_OF_ PHOSPHATIDYLINOSITOL, 3_KINASE SIGNA 
LING, KEGG_ERBB_SIGNALING_ PATHWAY, 
GO_ENDOTHELIAL CELL_ PROLIFERATION 
GO_CYSTEINE_TYPE_PEPTIDASE_ ACTIVITY GO_ENDOPEPTIDASE_ACTIVITY, 
GO_OXIDOREDUCTASE _ ACTIVITY_ACTING_ON_PAIRED_DONORS_WITH_INC 
ORPORATION_OR_REDUCTION_OF MOLECULAR_OXYGEN 
GO_DNA_BINDING TRANSCRIPTION_ACTIVATOR_ACTIVITY, 
REACTOME_NEGATIVE_REGULATION_OF_THE_PI3K_AKT_NETWORK, 
GO_POSITIVE_REGULATION_OF_ PHOSPHATIDYLINOSITOL, 3_KINASE SIGNA 
LING 
KEGG_ERBB_SIGNALING PATHWAY, MODULE_197, 
GO_COFACTOR_METABOLIC PROCESS 


KEGG_ERBB_SIGNALING_ PATHWAY, MODULE_197, 
GO_COFACTOR_METABOLIC PROCESS 


E"D" 表 示 该 化 合 物 为 初始 输入 药物 
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